High-energy gamma-ray afterglows from low-luminosity gamma-ray bursts 
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ABSTRACT 

The observations of gamma-ray bursts (GRBs) such as 980425, 031203 and 
060218, with luminosities much lower than those of other classic bursts, lead to the 
definition of a new class of GRBs - low-luminosity GRBs. The nature of the out- 
flow responsible for them is not clear yet. Two scenarios have been suggested: one 
is the conventional relativistic outflow with initial Lorentz factor of order of To > 10 
and the other is a trans -relativistic outflow with Fq — 1-2. Here we compare the 
high energy gamma-ray afterglow emission from these two different models, taking 
into account both synchrotron self inverse-Compton scattering (SSC) and the external 
inverse-Compton scattering due to photons from the cooling supernova or hypemova 
envelope (SNIC). We find that the conventional relativistic outflow model predicts a 
relatively high gamma-ray flux from SSC at early times (< lO'^s for typical parame- 
ters) with a rapidly decaying light curve, while in the trans-relativistic outflow model, 
one would expect a much flatter light curve of high-energy gamma-ray emission at 
early times, which could be dominated by both the SSC emission and SNIC emission, 
depending on the properties of the underlying supernova and the shock parameter eg 
and tB- The Fermi Gamma-ray Space Telescope should be able to distinguish between 
the two models in the future. 

Subject headings: gamma-rays: burst 



1. INTRODUCTION 

Long duration gamma-ray bursts are generally believed to result from the death of massive 
stars, and their association with core-collapse supemovae (SNe of Type Ib/c) has been observed 
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over the last decade. The first hint for such a connection came with the discovery of a nearby SN 
1998bw in the error circle of GRB 980425 (Galama et al. 1998; Iwamoto et al. 1998) at distance 
of only about 40Mpc. The isotropic gamma-ray energy release is of the order of only 10"^^erg 
(Galama et al. 1998) and the radio afterglow modelling suggests an energy of 10"*^- 10^"erg in a 
mildly relativistic ejecta (Kulkami et al. 1998; Chevalier & Li 1999). Recently Swift discovered 
GRB 060218, which is the second nearest GRB identified so far (Campana et al. 2006; Cusumano 
et al. 2006; Mirabal & Halpem 2006; Sakamoto 2006). It is also an under-energetic burst with 
energy in prompt 7/X-rays about 5 x 10"^^ erg and is associated with SN 2006aj. Another burst, 
GRB 031203, is a third example of this group (Sazonov et al. 2004). 

If one assumes that GRB 060218-like bursts follow the logN-logP relationship of high-luminosity 
GRBs then, as argued by Guetta et al. (2004), no such burst with redshifts z < 0.17 should be ob- 
served by a HETE-like instrument within the next 20 years. Therefore, the unexpected discovery of 
GRB 060218 may suggest that these objects form a different new class of GRBs from the conven- 
tional high-luminosity GRBs (Soderberg et al. 2006; Plan et al. 2006; Cobb et al. 2006; Liang et 
al. 2006; Guetta & Delia Valle 2007; Dai 2008), although what distinguishes such low-luminosity 
GRBs (LL-GRBs) from the conventional high-luminosity GRBs (HL-GRBs) remains unknown. 
Their rate of occurrence may be one order of magnitude higher than that of the typical ones (e.g. 
Soderberg et al. 2008). They are rarely recorded because such intrinsically dim GRBs can only 
be detected from relatively short distances with present gamma-ray instruments. 

The "compactness problem" of GRBs requires that the outflow of normal GRBs should be 
highly relativistic with a bulk Lorentz factor of Fq > 100 (Baring & Harding 1997; Lithwick & Sari 
2001). The low-luminosity and softer spectra of low-luminosity GRBs relax this constraint on the 
bulk Lorentz factor. It was suggested that the softer spectrum and low energetics of GRB060218 
(or classified as X-ray flash 060218) may indicate a somewhat lower Lorentz factor of the order of 
~ 10 (e.g. MazzaU 2006; Fan et al. 2006; Toma et al. 2006). 

An alternative possibility is that low-luminosity GRBs are driven by a trans-relativistic out- 
flow with To ~ 2 (Waxman 2004; Waxman, Meszaros & Campana, 2007; Wang, Li, Waxman & 
Meszaros 2007; Ando & Meszaros 2008). The flat light curve of the X-ray afterglow of the nearest 
GRB, GRB 980425/SN 1998bw, up to 100 days after the burst, has been argued to result from the 
coasting phase of a mildly relativistic shell with an energy of a few times 10"^^ erg (Waxman 2004). 
From the thermal energy density in the prompt emission of GRB060218/SN2006aj, Campana et 
al. (2006) inferred that the shell driving the radiation-dominated shock in GRB 060218/SN 2006aj 
must be mildly relativistic. This trans-relativistic shock could be driven by the outermost parts of 
the envelope that get accelerated to a mildly relativistic velocity when the supernova shock accel- 
erates in the density gradient of the envelope of the supernova progenitor (Colgate 1974; Matzner 
& McKee 1999; Tan et al. 2001), or it could be due to a choked relativistic jet propagating through 
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the progenitor (Wang et al. 2007). 

In this paper, we investigate the high energy afterglow emission from low luminosity GRBs 
for both relativistic and trans-relativistic models and explore whether the Fermi Gamma-ray Space 
Telescope can distinguish between these two models with future observations of low-luminosity 
GRBs. Since photons from the underlying supernova are an important seed photon source for 
inverse-Compton (IC) scattering, we consider both synchrotron self-inverse Compton (SSC) and 
external IC scattering due to supernova photons (denoted by SNIC hereafter). At early times 
(within a few days after the burst), a UV-optical SN component was recently detected from SN2006aj 
and SN2008D (Campana et al. 2006; Soderberg et al. 2008), which has been interpreted as the 
cooling SN envelope emission after being heated by the radiation-dominated shocly (Waxman et 
al. 2007; Soderberg et al. 2008; Chevalier & Fransson 2008). In our calculation, we take into 
account this seed photon source in addition to the late-time supernova emission, which peaks after 
ten days. 

Recently, Ando & Meszaros (2008) discussed the broadband emission from SSC and SNIC 
for a trans-relativistic ejecta in a low-luminosity GRB at a particular time- the ejecta deceleration 
time. Here we study the time evolution of the high-energy gamma-ray emission resulted from 
such IC processes and consider both the trans-relativistic ejecta and the highly relativistic ejecta 
scenarios for low-luminosity GRBs. 

The paper is organized as follows. First, we describe the dynamics of shock evolution in 
the two models in § [2l In §3, we present the formula for the calculation of the inverse-Compton 
emissivity. For the SNIC emission, we take into account the anisotropic scattering effect. Then we 
present the results of the spectra and light curves of SSC and SNIC emission for the two different 
models and explore the detectability of these components by Fermi Large Area Telescope (LAT) 
in § 4. Finally, we give the conclusions and discussions. 



2. Dynamics and Electron Energy Distribution 

We assume that the two models have the same parameters except for the initial Lorentz factor 
of the ejecta. For the latter, we adopt nominal values of To = 10 in the conventional relativistic 
ejecta model and Fq = 2 in the trans-relativistic ejecta model, respectively. Note that in the conven- 
tional relativistic ejecta model, even if Fq S> 10 the dynamics of the blast wave is identical to the 
case of Fq = 10 from the time tens of seconds after the burst, because the blast wave has entered 



'Although there is a disagreement on the origin of the early UV-optical emission from SN2006aj, an agreement 
has been reached for that of SN2008D (Waxman et al. 2007; Soderberg et al. 2008; Chevalier & Fransson 2008). 
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the Blandford-McKee self-similar phase since then. 

We consider a spherical GRB ejecta carrying a total energy ofE = lO^^E'soerg expanding into a 
surrounding wind medium with density profile n = Kr~^, where K = M/{A'KmpV^) = 3 x 10^^cm"^m 
with m = (M/10"^A/Qyr"^)/(Vw/10^kms"'). As the circumburst medium is swept up by the blast 
wave, the total kinetic energy of the fireball is (Panaitescu et al. 1998) 

£k = (7-l)("^ej+m)c' + 7C/' (1) 

where 7 is the bulk Lorentz factor of the shell, mej the ejecta mass, m is the mass of the swept- 
up medium, and the comoving internal energy U' can be expressed by i7' = (7- l)mc^, which is 
suitable for both ultrarelativistic and Newtonian shocks (Huang et al. 1999). Hereafter superscript 
prime represents that the quantities are measured in the comoving frame of the shell. As usual, we 
assume that the magnetic field and the electrons have a fraction eg and eg of the internal energy, 
respectively. Following Eq. ([T]), the differential dynamic equation can be derived as (Huang et al. 
2000) 

^ = - ^'-^ (2) 
dm niej + 27m ' 

which describes the overall evolution of the shell from relativistic phase to non-relativistic phase. 
The initial value of 7 is Tq = E /(niejC^). To obtain the time-dependence of 7 one makes use of 



dm 



= AnR^nm^, (3) 



dR 

— = /?c7(7+V7^), (4) 
dt 

where t is the observer time, R and j3 = 7"^ are the radius and velocity of the shell, respec- 
tively. 

Solving Eqs. 2, 3, and 4, three dynamic phases can be found: (i) Coasting phase. The shell 
does not decelerate significantly until it arrives at the deceleration radius, Rdcc = E / {A-KKV^m^c^), 
where the mass of the swept-up medium m is comparable to mej/Fo (Sari & Piran 1995). The 
corresponding deceleration time can be calculated from ?dec — ^dec/(2rQc). For representative 
parameter values £50 = 1 and m = 1, 7?dec — 4.4 x lO^^cm and ?dec — 1.8 x lO'^s for a trans-relativistic 
ejecta with Fq = 2, while 7?dec — 1.8 x lO'^^cm and ?dec — 29s for a conventional highly relativistic 
ejecta with Fq = 10; (ii) Blandford & McKee self-similar phase {t > fdec and 7 > 2), where 7 oc r^^^ 
and R oc t^^^; (iii) Non-relativistic phase (7 — > 1), where /? oc r^^^ and R oc t^^^, i.e., the Sedov-von 
Neumann-Taylor solution applies (Zel'dovich & Raizer 2002, p. 93, Waxman, 2004). 

In the absence of radiation losses, the energy distribution of shock- accelerated electrons be- 
hind the shock is usually assumed to be a power-law as dNg/d-ye oc As the electrons are cooled 
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by synchrotron and IC radiation, the electron energy distribution becomes a broken power-law, 
given by 



(1) for7g,c < le 



(2) for 7^,n, < 7^,c < le 
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7, 
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le.m <le< 7e,c 
77"^ 7.,c < 7.. < 7^;,n 



(6) 



which are normalized by the total number of the electrons solved from the dynamic equations. The 
minimum, cooling, and maximum Lorentz factors of electrons are, respectively, given by 



7e,m 
7e,c 



'p-1 m. 



(7-l) = 92/„ie._o.5(7-l) 
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(9) 
(10) 



where fpi = 6{p-2)/{p- 1), B' is the comoving magnetic field strength and Y is the Compton 
parameter that is defined as the ratio of the IC luminosity (including SSC and SNIC) to the syn- 
chrotron luminosit>@. 



3. Inverse-Compton Emission 

The accelerated electrons can be cooled by synchrotron radiation and inverse Compton (IC) 
scattering of seed photons (including synchrotron photons and blackbody photons emitted by the 
supernova). Since the IC emissivity on the basis of the Thomson cross section is inaccurate for 
high-energy 7-rays, we use the full Klein-Nishina cross section instead. Once the electron distri- 
bution and the flux of seed photons (/^,) (the distribution of which is isotropic) are known, the IC 



^At very late times, when ^^^^ decreases to be close to a few, ^e,m = ee^^(7~ 1)+ 1 is used (Huang & Cheng 
2003). 
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emissivity (at frequency v') of electrons can be calculated by(Blumenthal & Gould 1970; Yu et al 
2007) 

-'IC/,,/ 



- — / az/j. ^ 



£ iso('^ ) = 3^7 / dje— I du- ^^(x,};), (11) 



where 7e,mi„ = max[7e,c,7e,m,^i^V("^eC^)], ^i,mm = 1^' m^c^ / Alj^i-f^ifi^c^ - hu')], X = 47e/^z/,7meC^ 3; = 
hu' /[xi'jemeC^-hi'')], and 

g(x,y) = 2y\ny + (l + 2y)(l-y) + l-^^(\-y). (12) 

2 (1+Jcy) 

In our case, because the radius of the supernova photosphere is much smaller than that of the 
GRB shock, supernova seed photons can be regarded as a point photon source locating at the center 
in the comoving frame of the GRB shock. These soft photons from supernova impinge on the shock 
region basically along the radial direction in the rest frame of the shock, so the scatterings between 
these photons and the isotropically-distributed electrons in the shock are anisotropic. For a photon 
beam penetrating into the shock region where the electrons are moving isotropically, the inverse 
Compton scattering emissivity of the radiation scattered at an angle 9sc relative to the direction of 
the photon beam in the shock comoving frame is (Aharonian & Atoyan 1981, Brunetti 2000, Fan 
et. al. 2008): 

£'AiC(z/',cos^sc) = ^ r'^'^dj/-^ 

J<m.„ -rl'^s' 2(1-0 boa-0 blu-o^^ 

where ^ = hh'' /{%meC^), bo = 2(1 -cos^sc)7e^'^s/("^e<^^) ^i^d hu!. <^ hu' <C %meC^bg/(l+b0). On 
integration over 9sc for whole solid angle (i.e. in the case that the photon distribution is also 
isotropic), Eq.(fT3l) reduces to Eq.(fTT]). i.e. the usual isotropic inverse Compton scattering emis- 
sivity. In the observer frame, the angle 6 between the injecting photons and scattered photons 
relates with the angle ^sc in the comoving frame by cos6'sc = (cos^-/3)/(l-/3cos6'), where P is the 
velocity of the GRB shock. 

The observed SSC and SNIC flux densities at a frequency u are given respectively by (e.g. 
Huang et al. 2000, Yu et al. 2007) 

^ssc^ r ^_§M^D'sm9de, (14) 
" Jo ^T^Dl 

and 

F-=^'f:^5!^^^tf.i„«, (15) 

where Dl is the luminosity distance and D = [7(1 -/5cos6')]"^ is the Doppler factor. Note that 
cos6'sc in Eq.(15) can be transformed to cos^ through cos6'sc = (cos^-/3)/(l -/?cos6'). 
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3.1. Seed Photons from the Supernova 

For seed photons from the supernova, we consider the contributions from two components. 
One is the early thermal UV-optical emission from the cooling supernova envelope after being 
heated by the radiation-dominated shock (Waxman et al. 2007; Chevalier & Fransson 2008). Such 
UV-optical emission has been observed recently by Swift UVOT from SN2006aj (Campana et al. 
2006) and SN2008D (Soderberg et al. 2008). Another component is the supernova optical emission 
at later time, powered by the radioactive elements synthesized in supemovae. 

The characterization of the emission from the cooling supernova envelope is expressed as 
follows, which is correct at ? > lO^s (Waxman, Meszaros, & Campana 2007): 

^0.02 

^--^^^^^• ^mJm:)0.03 o^"^^- (1^) 

Lph(0 = 47ri?ph(0Vrph(0', (18) 

where i?ph(0 and T^hit) are the radius and temperature of the envelope, Lph(0 is the luminosity from 
the photosphere of the cooling envelope, Msnej and E^^^^ are the supernova ejecta mass and energy, 
and Rq is the initial stellar radius of the SN progenitor. We take the following values for supemovae 
like SN2006aj and SN2008D: £'snej,5i = £:snej/(10^^erg) = 2, Ms„ej = IMq and i?o,i2 = i?o/(10^^cm) = 
0.3 (Mazzali et al. 2006; Soderberg et al. 2008). For hypemovae such as SN1998bw associated 
with GRB980425, modelling of the SN optical emission gives a larger kinetic energy and ejecta 
mass such as £'s„ej,5i ~ 22 and Ms„ej = 6M0 for SN1998bw (Woosley et al. 1999). 

The luminosity of the late component is rising before ~ lOdays and then decaying exponen- 
tially. For SN2008D, it is found that the rising is roughly in proportion to (Soderberg et. al. 
2008), so we can describe the luminosity as 



_. 3 X 10'*2(^^)i-6erg s-\ ? < lOdays, 

W0-<; 3xl0^2,^;(^;__^)^^g^-i^ r>10days ^^^^ 



We assumed that the radiation temperature is approximately constant, r^d — leV. The bolometric 
luminosity of the cooling supernova envelop is dominant over that of the late component before 
f = 7.8 X lOV 

Through the calculation of the shock dynamic evolution, the radii of the GRB shock are 
/? = 1.9 X 10^"^ cm at ? = lO^s, and i? = 1.5 x lO^^cm at ? = lO'^s for the trans-relativistic ejecta 
model, and for the conventional ejecta model = 1.0 x lO^^cm at ? = lO^s and R = A. \ x lO^^cm 
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at t = lO^s. According to Eq{T6l the shell radii for supernova are 9.6 x lO^^cm at t = lO^s and 
6.1 X lO'^cm at f = lO'^s, and the one for hypernova are 1.8 x lO^^cm att = lO^s and 1.1 x lO'^^cm 
att= lO'^s. Consequently the shell radius is less than ten percents of the GRB ejecta height, so it 
is reasonable to approximate that supernova photons come from a point source at the center and 
impinge onto the electrons in the shock from behind. 



3.2. Compton Parameter Y 

Following Moderski, Sikora & Bulik (2000) and Sari & Esin (2001), we define the Compton 
parameter Y as 

y _ Lie _ Mj^ _ ^^SYN + /a^^SN 

^SYN Mb ^ b 



where 



/ -2 -^'SN 

MSN = 7 



65 erg cm-X^ %jr\ (21) 



M SYN 



IjeeU' 



(1+7) 

^ 5 . 4 X 1 0" erg cm"^ e,,_,,5 mRlj^i-^ -IW', (22) 
u'b = esu' 

~ 1.8 erg cm-\B,_3mi?727(7- 1), (23) 

are, respectively, the comoving energy densities of the blackbody supernova seed photons, syn- 
chrotron seed photons and magnetic fields, and fa is the factor accounting for the suppression of 
the photon energy density due to the anisotropic inverse Compton scattering effect (/a = 1 corre- 
sponds to the isotropic scattering case). Here Lsn = ^ph + ^rad^ which is dominated by the lumi- 
nosity of the cooling envelop Lph, u' is the comoving internal energy density, rj = r/rad^KN is the 
radiation efficiency where rj^^d is the fraction that the electron's energy radiated, and t^kn is the 
fraction of synchrotron photons below the KN limit frequency (Nakar 2007). For slow cooling, 
r?rad = (le,c/le,mf~'' (Modcrskl, Sikora & Bulik 2000), and 

i^^cn(7..c) < < 




VKN={ (^^)W, U'^<ui^hej<u'^ (24) 

K < i^U(7«,c) 



For fast cooling, rj^-ad = 1 and 



Vkn 



0, ^kN(le.m) < l^c 

c-H7^y^\ K<^Ui.^)<< (25) 

1, < < I^^CN(7.,m) 
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Solving eq. (I20l) . we get 



y = i 

2 



+ 



SN 



-1 



Roughly, the above expression can be simplified in three limiting cases as follows 



Y = 



(A/477e,/es+l- l)/2, /aW'sN < min[M'sYN,MB] 



(26) 



(27) 



In the latter two cases, Y can be treated as a constant when the electrons are in the fast-cooling 
regime and u'^ < u^i-fe.m) with 77 = 1. 



3.3. Pair Production Opacity For High Energy Photons 

High-energy gamma-rays can be attenuated due to interaction with low-energy photons through 
the pair production effect. We consider the pair production opacity in the shock frame due to the ab- 
sorption by low-energy photons, which include thermal photons from the supernova, synchrotron 
photons, SSC photons and SNIC photons. A high energy photon of energy E'^ j in the shock frame 
will annihilate with a low energy photon of provided that E'^ ^E^2i^~^^^^n) > 2(meC^)^, 
where is the collision angle of the two annihilation photons. The pair creation cross section is 
given by 

a(£;i,£;2) = ^7rr2(l-^2~,^(3_^4~)^j^ll|_2^(2-^2-,j, (28) 

where (3 = v/c= ^Jl -2(mgC^)2/[£':^ j£'^2(l~'^os^i2)] is the velocity of electrons in the center-of- 
mass frame (Heitler 1954, Stecker, De Jager, & Salamon, 1992) and ro = /{mgC^) is the classic 
electron radius. By using the photon distribution in the shock frame, we obtain the optical depth 
for high-energy photons of energy, 

T(E'^,) = / aiE'^„E'^2h-yiE'^2)-dE'^2, (29) 

where r]R is the shock compressed ratio, which is ?7r = 47-1-3, and the threshold energy of low 
energy photons is E^^^ = lim^c^)^ / E'^ j(l -cos 6*12). In the calculation, we estimate the cutoff energy 
conservatively by assuming that colliding photons are moving isotropically in the shock frame. 
Such a treatment may overestimate the pair-production opacity because in reality the supernova 
seed photons move anisotropically (i.e. moving outward in radial direction as seen by high-energy 
photons emitted from the shock at much larger radii). 
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The result of the cutoff energy in the observer frame fi^.cut = lE'^if = 1) is given in Figure [B 
Since the luminosity and peak energy of the supernova envelope emission decreases with time in 
general and the shock radius increases with time, the cutoff energy of the high energy spectrum 
increases with time, which is clearly seen in Fig. [T] From this figure, we can see that the cutoff 
energy is above IGeV after the starting time of our calculation (lO^^s) , so we can calculate the 
light curves at energy ~ IGeV without considering the opacity. Since the cutoff energy is large 
enough, it does not affect the detectability of Fermi LAT, whose sensitive energy band is 20MeV ~ 
300GeV. 



4. Results 

4.1. Spectra and Light Curves: Numerical Results 

We first compare the spectra of different IC components for the two models. The spectra at 
t = 10'' s are shown in Figured including the synchrotron emission, the thermal emission from 
the supernova or hypernova, the synchrotron self-Compton emission and the SNIC emission, for 
parameters e^, = 0.1, = 0.001, p = 2.2, E = 10^°erg and burst distance Dl = lOOMpc. From 
the spectra we can see that for the supernova case the SSC emission dominates over the SNIC 
emission at energies from IMeV to lOOGeV for conventional relativistic ejecta model, and for 
trans-relativistic ejecta model the SSC emission is higher than SNIC emission and the two com- 
ponents are both important at high energies below the cutoff energy. For the hypernova case the 
SSC emission is dominant in the conventional relativistic ejecta model while the SNIC emission 
is dominant in the trans-relativistic ejecta model at the high energy band. This is because in the 
case of hypemovae and Fq = 2, the energy density of hupemova photons (multiplied by the 
suppressed factor /a due to the anisotropic scattering) is significantly higher than the synchrotron 
radiation density Msyn- 

Figure [3] shows light curves of SSC and SNIC afterglow emission at hu = IGeV in the two 
models for E = 10^°ergs, p = 2.2 and four different sets of parameters of and eg for the supernova 
case. In the conventional relativistic model with Fq = 10, a sharp decay phase of a GeV afterglow 
is produced during the early hours, which is mildly dominated by the SSC emission, and a slightly 
flatter decay phase dominated by SNIC emission takes over at late times t > lO^s. On the other 
hand, for the trans-relativistic ejecta model with Fq = 2, a plateau of SSC emission, due to the 
presence of a coasting phase in the ejecta dynamic, dominates in the early time, which transits to a 
faster decay at later time and SNIC emission become dominant after the time 10^ ~ lO^s. 

Figure|4]show light curves for the hypernovae case. Light curves in the conventional relativis- 
tic ejecta model are similar with those for the supernova case except that SNIC emission become 
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dominant at earlier time; for the trans-relativistic ejecta model, SNIC emission is always domi- 
nated with parameters e^. = 0. 1 and eg = 0.001, a plateau is seen in the early time, which transits to 
a faster decay at later time. 

The two models also predict different flux levels at high-energies. In early hours, the total flux 
from the conventional relativistic ejecta model is more than one order of magnitude higher than 
that in the trans-relativistic ejecta model. This can be explained by the different amount of energy 
in shocked electrons. For the conventional relativistic ejecta model, the ejecta has been decelerated 
at this time and a great part of its energy has been converted into shocked electrons, while in the 
trans-relativistic ejecta model, only a small fraction of the ejecta kinetic energy has been converted 
to shocked electrons at this early time. A lower amount of energy in shocked electrons results in a 
lower flux level in the trans-relativistic ejecta model. 

In addition, a comparison among four panels in Figure [3] and Figure |4] indicates that the flux 

decreases as decreases, which is obvious since the energy of radiating electrons Ee oc e^. As eg 

decreases, the SSC flux changes little and the SNIC flux increases for the case that SSC emission 

is dominant at early time. This can be understood from the following analysis. Since the SSC 

-- - 

emission dominates and /aMsN ^ the early times, F ~ y rye^/ eg cx e^^ e| . The SSC and SNIC 
flux at z/ > z/min > i^c scalc as 

'^^f = ^^ma'x^(^)"K-ic)"^ « (30) 
c min 



and 



SNIC 1 3 

Ur^ -i^-^max V SNIC'' \ .SNIC ' OC 6^ Ce . pi; 



4.2. The effect of the anisotropic scattering on the SNIC emission 

The incoming supernova photons are anisotropic as seen by the isotropically distributed elec- 
trons in the GRB shock, so the IC scatterings are anisotropic. In order to see how the anisotropic 
inverse-Compton scattering (AIC) affects the SNIC flux, we compare the light curves of the SNIC 
emission obtained by using the isotropic scattering formula Eq. 1 1 and using the AIC scattering 
formula Eq.l3 in Figure [51 The thinner lines show the SNIC light curves obtained using the usual 
isotropic scattering formula, while the thicker ones correspond to the calculations with the AIC 
scattering effect taken into account. One can see that the flux of the SNIC emission with the AIC 
effect correction is reduced by a factor of about ~ 0.4 compared to the isotropic scattering case. 
This is consistent with the calculation result obtained by Fan & Piran (2006), who studied the 
anisotropic inverse Compton scattering between inner optical/X-ray flare photons and electrons in 
the outer GRB forward shock. 



-12- 



Fig.5 shows that the AIC effect suppresses the SNIC flux only slightly. The anisotropic photon 
distribution results in more head-on scatterings, i.e. the photon beam scatter preferentially with 
those electrons that move in the direction antiparallel to the photon beam, so one can expect that 
the scattered IC emission power has a maximum at 9sc = tt and goes to zero for small scattering 
angles (e.g. Brunetti 2000). The photons scattered into the angles < 9sc ^ 7r/2 relative to the 
shock moving direction in the shock comoving frame will fall into the cone of angle l/F in the 
observer frame, according to the transformation formula cos^sc = (cos9-/3) /(I -/3cos^). Therefore 
the AIC scatterings decrease the IC emission in the 1 /F cone along the direction of the photon 
beam, but meanwhile they enhance the emission at larger angles (about half of the emission falling 
into angles between l/F and 2/F, see Wang & Meszaros 2006). For a spherical outflow as we 
consider here, the IC emission after integration over angles should have the same flux in every 
direction in the observer frame, with a flux level only slightly reduced comparable to the isotropic 
scattering case. 



4.3. Analytical Light Curves 

As a comparison, we derive here approximate analytical expressions for afterglow light curves, 
which provide an explanation for the physical origin of the behavior. Since the anisotropic SNIC 
emission are depressed by a factor of 0.4, which is almost constant, relative to the isotropic seed 
photons case, we can consider the isotropic seed photons case for the approximate analytic treat- 
ment of the afterglow light curves. The blackbody photons from the supernova can be approxi- 
mated as mono-energetic photons with /i^sn = 2.7KTsn. Thus, similar to the description for syn- 
chrotron emission of a single electron (Sari et al. 1998), the radiation power and characteristic 
frequencies of SNIC from a single electron scattering supernova photons in the observer frame can 
be described by 

P(7.)=^W7'^^, (32) 

and 

u(-fe) = I'J'J^I'sn/I = (33) 



respectively. Similar to the analysis in Sari and Esin (2001), the maximum flux of the SNIC 
spectrum is 

^sNic^Jv^:^ .34. 

Characteristic SNIC frequencies are 

^iT = 2x',m^SN (35) 
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and 

^r'' = 27jc^sN (36) 

respectively. 

By adopting the broken power-law approximation for the IC spectral component ( Sari & Esin 
2001) and the dynamics of the shock discussed in § |2l one can derive the analytic light curves in 
an approximate way. 



4.3.1. Light Curves In The Conventional Relativistic Ejecta Model 
In the conventional relativistic ejecta model, at time ?dec < t < ti^^a, where ?tran is defined as 



the transition time when 7^ m = 7^ c, we have 7]=! and u'^ <C /a^sN ^ "syn' so 7 ~ A/r/e^/eB oc 
The shock dynamic follows the Blandford & McKee self-similar solution in the wind medium, i.e. 
i? oc ?2 and 7 oc Then we can obtain the evolution of the break frequencies of the SSC and 
SNIC spectral components and their peak flux in the following way: 

^Hn = OC l^f ^ = Ij^U, OC t\ (37) 

pSSC _ ^tNc sYN .-1 /oox 

and 

^r"«r\z.r"«^F,^,rx«r«-'- (39) 

where z/sn o<i Tsn oc t~^/'^ has been used for the cooling envelope emission (see Eq. 13). The SSC and 
SNIC flux at an observed frequency v higher than characteristic frequencies vary as vF^^'^ oc t~''*'^ 
and vF^^^*^ oc rz'^" *' respectively for t < ?tran- 

At f > fsNic, where fsNic is the time when f^u'^^ = u'^yj^, we take Y ~ /aw'sN/(eB"') because the 
SNIC emission becomes dominated. At such time, the shock is likely to enter the non-relativistic 
phase, so we take the Sedov-von Neumann-Taylor solutions 7? oc fa and (3 (xt~^, which induce that 
Y oc ? 3"° "^. So the minimum and cooling Lorentz factors of electrons vary as 7e m oc oc Ta and 
7e c OC (3~^R^t~^Y~^ oc t'^'^*^. Thus, the break frequencies of the SSC and SNIC spectral components 
and their peak fluxes evolve with time in the following way: 

u'J^ oc TT, zyf c ^ oc r\ (40) 



and 



ocr^,^r' ocr^-o^Fj;:^, ocrt-°^ (41) 



Then the SSC and SNIC fluxes at high energy u vary as vF^^^ oc '« ^ and vF^^^^ oc ?6 '2'' at 

^ > ^SNIC- 



-14- 



To summarize, the temporal evolution of the SSC and SNIC afterglow emission at high ener- 
gies are 

zvFf c oc I 2 ^ ^ ^'^^ (42) 



and 

in the two asymptotic phases for p = 2.2. 



4.3.2. Light Curves In The Trans-relativistic Ejecta Model 

In the trans-relativistic ejecta model, one would expect a much flatter light curve of high- 
energy gamma-ray emission at early times, which could be dominated by both the SSC emission 
and SNIC emission, depending on the properties of the underlying supernova and the shock pa- 
rameter te and eg. For the supernova case, the SSC emission is dominant before the deceleration 
time and the transition time for most parameters, while for the hypemova case, the SSC emission 
is dominant at earlier time with = 0.3 and eg = 0.01 and the SNIC emission is always dominant 
for 6^ = 0.1 and eg = 0.001. 

For cases where SSC emission dominated in trans-relativistic ejecta model, we have = 1 in 
fast cooling regime and Mg <^ /a^sN ^ "syn' so 7 ~ \fr\^7[^B t^. Since t < ?dec> we adopt the 
approximation R(xt,^(xt°. So the SSC and SNIC flux vary as uF^^'^ oc t^~2 and uF^^'^ oc t°-'^~^ 
respectively att < min(?dec5 ^tran)- At later time when t > t^ec and 7 1, the SNIC emission become 
dominant, the evolution of light curves is the same as that in conventional relativistic ejecta model, 
i.e. the SSC and SNIC flux vary as uF^^'^ oc t^-^-TP and i^F^^^'^ oc f-^'riP at t > ^dec- Therefore, 
for the supernova case and the hypernova case with early dominated SSC emission, the temporal 
evolution of the SSC and SNIC emission at high frequency uforp- 2.2 are given by 



^.Ff^oc<' '_2.2 (44) 



l ? <min(?dec,?tran) 

_1 O 



0C< 

\^ t t ^ tdec 

and 



j.TpSmC ^ \ ^^'^^ ^ ^ niin(?dec)^tran) 

oci ^_i.2 ^^^^^ (4:)) 



For cases where SNIC emission dominated in trans-relativistic ejecta model. We have /aw'sN ^ 
m'syn' so y ~ fafi'sti/i^Bu'). At the time t < ?dec> we adopt the approximation R (xt, ^ oat^ Xo yield 
Y oc f^-^. So the SSC and SNIC flux vary as vF^^^ oc t^-^'^ and vF^^^^ oc respectively at 
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t < min(?ciec5^tran)- At later time when t > tdec and 7 — > 1, the evolution of light curves is the same as 
that in conventional relativistic ejecta model, i.e. the SSC and SNIC flux vary as uF^^^ oc t^-^~~6P 
and z/Fjf^^^ oc t^'nP att > ?dec- Therefore, for the supernova case and the hypemova case with early 
dominated SSC emission, the temporal evolution of the SSC and SNIC emission at high frequency 
for p = 2.2 are given by 

ZvFf C oc I ''22 ' ~ l^in^^dec, ^tran) ^^^^ 



and 



^.Ff^^oc<! ^-"^^^^^'i-'^^-^ (47) 



t ' t ^ ?dec 



4.4. Detectability by the Fermi LAT 

We explore here whether Fermi LAT can detect the high energy gamma-ray emission from low 
luminosity GRBs in the two models considered above. Following Zhang & Meszaros (2001), Gou 
& Meszaros (2007) and Yu, Liu & Dai (2007), the fluence threshold for long-duration observations 
is Fthr = [(t)o(t I ?eff)'''^]£^ph^eff wWch is in proportional to {"1^ due to the limitation by the background, 
where we take the average energy of the detected photons as E^<^ = 400MeV and the effective time 
as ten = 0.5yr. (po is the integral sensitivity above 100 MeV for LAT for a steady source after a year 
sky survey, which is 0o ~ 3 x 10"^phscm"^s"^ (atwood et. al. 2009) and is improved by a factor of 
3 by keeping the GRB position at the center of the LAT field of view as long as possible (Gou & 
Meszaros 2007). For short-time observation, the fluence threshold is calculated by Fthr = 5£'ph/Aeff 
under the assumption that at least 5 photons are collected. Taking the effective area Aeff = 6000cm^, 
we can obtain the fluence threshold of Fermi LAT, 

J 5.3 X 10-Vg cm-2, ?<4.4xl0^s, 
" \ 2.5 X lO-Y/^erg cm-\ ? > 4.4 x lO^s. ^ ^ 

With this fluence threshold, the detectability of high energy emission (with supernova seed photons 
luminosity given above and a total energy E = 10^° erg) by the Fermi LAT is shown in Figure 6. 
The time-integrated fluence shown in the plot is defined as an integration of the flux density (F^) 
over the Fermi LAT energy band [20MeV, 300GeV] and the time interval [0.5?, as used in Gou 
& Meszaros (2007), which is J^ ^^ £^ F^dudt. 

For Se = 0.3 and other representative parameter values, high-energy gamma-ray emission in 
the conventional relativistic ejecta model stays detectable up to ~ lO^s, while the high-energy 
gamma-ray emission in the trans-relativistic ejecta model can only be detected in a short period 
around lO'^^s. For a lower value such as = 0.1, the high-energy gamma-ray emission can still 
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be detected in the conventional relativistic ejecta model, while it becomes undetectable for the 
trans -relativistic ejecta model. 

In order to compare with earlier results in Ando & Meszaros (2008), we increase the total 
energy to £ = 2 x 10^°erg, which yields a kinetic energy E]^ = (Tq- l)/ToE = 10^°erg in the trans- 
relativistic ejecta model, the same as that used in Ando & Meszaros (2008). This will increase the 
fluence by a factor of 2. We also increase the SN luminosity from that of a normal SN (shown 
in § 13.11 ) to SN1998bw-like hypemovae. In Figure 7, we show the detectability of high energy 
emission by Fermi LAT in this case. By comparing the flux of light curves between the supernova 
case and the hypemova case, which is shown in Fig. |3]and Fig. HI we can see that increasing the 
SN luminosity can hardly enhance the total IC flux for the two models. Our IC flux is still lower 
than that obtained by Ando & Meszaros (2008). The main difference between our calculation and 
that of Ando & Meszaros (2008) is the different minimum Lorentz factors used in the calculations. 
Ando & Meszaros (2008) may overestimate the minimum Lorentz factor of electrons by taking 
7e,m = '^eimp/nie)^, wMch is a factor of {p-\)/{p-2) larger than ours. From the formula describing 
the high energy flux vF^^^'^ oc t'mta snic ^ 7e'm ' oi^^ expects that the flux is increased by a factor of 
which is about 8 for p = 2.2. 



5. Discussions and Conclusions 

The external stellar wind provides a source of Thomson opacity to scatter the supernova emis- 
sion, thus a quasi-isotropic, back- scattered SN radiation field is present. Let's study whether this 
component is important. For a GRB shock locating at radius R and moving with a Lorentz factor 
7, the Thomson scattering optical depth of the wind is = ajRn = ajKR'^ = 2.0 x lO~'^mRll. 
The scattered SN energy density by the wind in the comoving frame of the blast wave is u'^ = 
{LsN /4-'ivcR^)T^'y^ due to the relativistic boosting effect. On the other hand, the energy density of 
the supernova photons impinging the shock from behind is MgN - I'^i^sN 1 47rc7?^) (see Eq.2 1 ). The 
ratio between these two energy densities is Msn/^sn = 2.0 x 10"^mi?757'^. For the trans-relativistic 
ejecta model with To = 2, the radius of the GRB ejecta is i? = 1 .9 x 10^^ cm at ? = lO^'s, so ^ "sn 
for typical wind parameters. For the conventional relativistic ejecta model, the Lorentz factor and 
radius of the GRB ejecta are respectively 7^4 and i? = 1.0 x 10^^ cm at f = lO^s, so we also have 
"sN ^ "sN- This means that the wind- scattered supernova radiation field is subdominant com- 
pared to the direct impinging supernova photon field and hence we neglect its contribution to the 
high-energy gamma-ray emission. 

Our estimate of the pair-production opacity for high-energy photons in § I3.3l is based on the 
common assumption that the colliding photons are isotropic in the rest frame. However, as we 
have shown above, the low-energy photons from the supernova essentially move radially outward 
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before colliding with high-energy photons (e.g. Wang, Li & Meszaros 2006). Therefore, the col- 
lision process between the high-energy photons and soft supernova photons is anisotropic, which 
would decrease the pair-production opacity. This will be the subject of a more detailed future cal- 
culation and would be useful to explore whether TeV photons can escape from the source, which 
is important for checking the detectability by ground-based Cherenkov detectors such as Magic, 
VERITAS, Milago, HESS, ARGO etc. Additionally, besides the high-energy gamma-ray emission 
discussed in this work, the high-energy neutrino emission arising from p'y interactions between 
shock- accelerated protons and photons from the supernova may also provide a constraint on the 
model for low-luminosity GRBs (e.g. Yu et al. 2008). 

Trans-relativistic ejecta may also exist in the usual high luminosity long GRBs, besides in low 
luminosity ones, since accelerating shocks are expected to accompany the supernova. Berger et al. 
(2003) and Sheth et al. (2003) found that the radio and optical afterglow indicates a low velocity 
component more than l.Sdays after the explosion in GRB030329/SN2003dh. However, since the 
highly relativistic ejecta is much more energetic than the trans-relativistic component, high-energy 
gamma-ray emission from the latter component could easily remain hidden. 

In summary, we have calculated the spectra and light curves of the high-energy gamma-ray 
afterglow emission from low luminosity GRBs for the two main models in the literature, i.e. 
the trans-relativistic ejecta model (Fq — 2) and the conventional highly relativistic ejecta model 
(Fo > 10), considering both synchrotron self inverse-Compton (SSC) and the external inverse- 
Compton due to photons from the underlying supemova/hypemova. Our analysis takes into ac- 
count a full Klein-Nishina cross section for inverse Compton scatterings, the anisotropic scattering 
of supernova photons and the opacity for high energy photons due to annihilation with low-energy 
photons (mainly from the supemova/hypemova). 

We find that for the supemova case the conventional relativistic outflow model predicts a 
relatively high gamma-ray flux from SSC at early times (< 10^ - lO^s for typical parameters) with 
a rapidly decreasing flux, while in the trans-relativistic outflow model, a much flatter light curve 
of high-energy gamma-ray emission dominated by the SSC emission, is expected at early times. 
For the hypernova case, the SSC emission also dominates and decays sharply at early time in 
the conventional relativistic ejecta model, while for the trans-relativistic ejecta model, both the 
SSC emission and the SNIC emission could be dominant at early time, depending on the shock 
parameters and eg. The main difference between these two models arises from their different 
initial Lorentz factors, which induces a different dynamical evolution of the shock at early times. 
As a result, different observational features arise, such as different light curve shapes, different 
flux levels and different dominant components (as detailed in § 14.11 ). As shown in § 14.41 high- 
energy gamma-ray emission can be detected in both models as long as is large enough, although 
detection from the conventional relativistic ejecta is much easier. Thus, with future high energy 
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gamma-ray observations by Fermi LAT, one can expect to be able to distinguish between the two 
models based on the above differences in observational features. 
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Fig. 1. — The high-energy cutoff due to annihilation with low-energy photons at different times in 
the two models for low-luminosity GRBs, with = 0.001, p = 2.2, m=l,E = lO^^erg and different 
parameters of e^. The solid lines and dashed lines show the cutoff energy in the conventional 
relativistic ejecta model with e^, = 0.3 and = 0.1, respectively. The dash-dotted and dotted lines 
show the cutoff energy in the trans-relativistic ejecta model with e^, = 0.3 and e^, = 0. 1, respectively. 
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Fig. 2. — The spectra of SNIC, SSC, synchrotron emission and supernova/hypernova photons in 
the two models at times t = lO^s for the parameters = 0. 1 , = 0.001 , p = 2.2 and E = 10^"erg. The 
left panels show the supernova case, the right ones show the hypemova case, the top ones denote 
the conventional relativistic ejecta model and the bottom ones denote the trans-relativistic ejecta 
model. The dotted and dash dotted lines show the spectra of SSC and SNIC with the annihilation 
effect taken into account, while the solid lines and dashed lines show the spcetra of SSC and SNIC 
emission without considering this effect. The short dashed lines denote the spectra of synchrotron 
emission and the short dotted lines denote the blackbody spectrum from the supernova/hypernova. 
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Fig. 3. — Light curves of the SNIC and SSC emission at energy hu = IGeV in the two models 
for a normal supernova case and for E = lO^^erg and different values for and eg. The thicker 
and thinner solid lines represent the SNIC and SSC emission in the trans-relativistic ejecta model, 
and the dashed and the dotted lines represent the SNIC and SSC emission in the conventional 
relativistic ejecta model. 
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Fig. 4. — Same as Fig. 3, but for the hypemova case. The thicker and thinner solid lines represent 
the SNIC and SSC emission in the trans-relativistic ejecta model, and the dashed and the dotted 
lines represent the SNIC and SSC emission in the conventional relativistic ejecta model. 
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Fig. 5. — Comparison of the light curves of the SNIC emission at hu = IGeV with (the thicker 
lines) and without (the thinner lines) the anisotropic scattering effect correction in the two models 
of low-luminosity GRBs. The solid lines represent the SNIC emission in the trans-relativistic 
ejecta model, and the dashed lines represent the SNIC emission in the conventional relativistic 
ejecta model. 
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Fig. 6. — Time-integrated fluence (defined as f^^^ F^dvdt, where hv^ = 20MeV and hu2 = 
min(i/cut:300GeV)) of the high-energy gamma-ray emission in the two models of low-luminosity 
GRBs. Short dotted lines represent the fluence threshold of Fermi LAT at 400MeV with effective 
detection area A^s = 6000cm^. The thicker and thinner solid lines represent the time-integrated 
fluence in the trans-relativistic ejecta model with = 0.3 and = 0.1, respectively; the dashed 
lines and dotted lines represent the time-integrated fluence in the conventional relativistic ejecta 
model with = 0.3 and = 0.1, respectively. Other parameteres are E = 10^°erg, m = 1 and 
eg = 0.001. 
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Fig. 7. — Same as Figl6l but for the case of a hypernova explosion that has a larger total energy 
(E = 2x lO^^erg), in accordance with the parameters used in Ando & Meszaros (2008). 



